q_hig=pars(1);
zlow=pars(2);
q_low=pars(3);
zhig=pars(4);
q_mid=pars(5);
zmid=pars(6);
rbar=pars(9);
r=pars(11);
delta=pars(13);
xi=pars(14);

coupon=pars(17);
rho=pars(18);
sigma=pars(19);
phi=pars(20);
a=pars(21);
varphi=pars(22);
nu=pars(23);
markupsub=pars(24);
m_ss=pars(26);
b_up=pars(27);
b_dn=pars(28);
b_md=pars(29);
theta=pars(32);
l_s=pars(33);

A=pars(34);
y=A^(1/(1-nu));
z_0=pars(12);
k=pars(7);
varzeta=pars(8);
tau=pars(10);
chi=pars(16);
B=pars(30);
R0=A/nu-chi;
zeta=pars(31);
zavg=pars(35);
%parameters that scale with value added
% A_cailb=pars(34);
% y_calib=pars(25);
% y=A^(1/(1-nu));
% z_0_calib=pars(12);
% z_0=z_0_calib*(y/y_calib)/(A/A_cailb);
% k=pars(7)*(A/A_cailb)*(z_0/z_0_calib);
% varzeta=pars(8)*(A/A_cailb)*(z_0/z_0_calib);
% tau=pars(10)*(A/A_cailb)*(z_0/z_0_calib);
% chi=pars(16)*(A/A_cailb);
% B=pars(30)*(A/A_cailb)*(z_0/z_0_calib);
% R0=A/nu-chi;
% zeta=pars(31)*(A_cailb/A)^(1/theta);
